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Abstract 

Effects of a generalized dark energy fluid is investigated on cosmic density fluc- 
tuations such as cosmic microwave background. As a general dark energy fluid, we 
take into consideration the possibility of the anisotropic stress for dark energy, which 
has not been discussed much in the literature. We comprehensively study its effects 
on the evolution of density fluctuations along with that of non-adiabatic pressure 
fluctuation of dark energy, then give constraints on such a generalized dark energy 
from current observations. We show that, though we cannot find any stringent limits 
on the anisotropic stress or the non-adiabatic pressure fluctuation themselves, the 
constrains on the equation of state of dark energy can be affected in some cases by 
the nature of dark energy fluctuation characterized by these properties. This may 
have important implications to the strategy to study the nature of dark energy. 



1 Introduction 



Almost all cosmological observations today suggest that the expansion of the present uni- 
verse is accelerating. To account for the present cosmic acceleration, one usually assumes 
an enigmatic component called dark energy. Although many authors have investigated the 
dark energy to date, we still do not know its nature yet. In most researches of dark energy, 
one often parameterizes dark energy with its equation of state wx = Px/px where px and 
px are pressure and energy density of dark energy. Current cosmological observations 
such as cosmic microwave background (CMB), type la supernovae (SNela), large scale 
structure (LSS) and so on can give a constraint on wx- When one puts constraint on the 
equation of state, it can be assumed to be constant in time for simplicity. However many 
models of dark energy proposed so far has a time-dependent equation of state. Thus most 
of recent study of dark energy accommodate its time dependence in some way. Although 
its time dependence can be complicated in general, it is usually assumed a simple form 
such as wx = wo + wi(l — a) [HE] when one tries to limit wx accommodating its time 
dependence. Some authors have also discussed the constraint on dark energy using the 
time evolution of the energy density of dark energy itself or Hubble parameter by binning 
to divide them into several epochs [3]. 

In addition to the time dependence of the equation of state, there is another property of 
dark energy of which we have to take care. In fact, in various models of dark energy, dark 
energy can fluctuate and then its fluctuation can affect the cosmic density perturbations 
such as CMB and LSS. When one specifies a model for dark energy, one can investigate 
its effects of fluctuation using the perturbation equations in the model. However, because 
there are many models proposed to date, one may prefer to describe it in a phenomenologi- 
cal way. To describe the nature of fluctuation of dark energy, the (effective) speed of sound 
c 2 s has been used in many works so far. Some authors have investigated the constraints 
on the equation of state varying the speed of sound [I] and on the speed of sound itself 
[SI El El- Since the fluctuation of dark energy mostly come into play after when dark energy 
becomes the dominant component of the universe, the effects of it appear on large scales. 
Hence, when one tries to limit the speed of sound from observations of CMB, one does not 
obtain a severe constraint on it because the cosmic variance error is large there. However, 
some studies have shown that the constrains on the equation of state can vary with the 
assumption for the perturbation property of dark energy. For example, when one assumes 
that there is no fluctuation of dark energy, the constraint for the constant equation of 
state is Wx = — Q-9^-oaoi I^J- However, when one takes into account the perturbation of 
dark energy, the constraint becomes wx = — l.OOj^g where c 2 s = 1 is assumed [8]. Thus, 
in this sense, it is important to investigate the properties of dark energy fluctuation from 
the viewpoint not only from to reveal its fluctuation nature but also from constraining the 
equation of state for dark energy. 

In fact, to parameterize the perturbation property of dark energy, the use of the speed 
of sound is not enough. If one considers the possibilities of imperfect fluid models for dark 
energy, one also needs to specify its anisotropic stress in some way. In fact, there are some 
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models proposed which have such anisotropic stress (for example, [9]). There are also a 
few works which accommodate such possibilities from a phenomenological point of view 
[TUl [TTl [T2] . In particular, the authors of Ref . [12] have discussed the implications of the 
anisotropic stress to observations such as CMB. Although they investigated its effects in 
some detail, cosmological constrains in such models was not given. Since we have now 
precise measurements of cosmology, it is worth to investigate constraints for dark energy 
in such a general setting. 

The purpose of the present paper is to discuss a generalized dark energy including 
the possibilities of anisotropic stress and give the constraints on cosmological parameters 
assuming such a generalized dark energy. 

The structure of this paper is as follows. In the next section, we briefly review the 
formalism for the investigation of a generalized dark energy. The perturbation equations 
are also given there. Then, in section 3, we discuss the effects of dark energy fluctuation 
on CMB, with particular attentions to the anisotropic stress of dark energy. In section 4, 
we present the constraints on dark energy from current cosmological observations. The 
final section is devoted to conclusions and discussions. 

2 Formalism 

In this section, we summarize the formalism to investigate a generalized dark energy. Some 
detailed descriptions of this issue can be found in Refs. [TQl [12] . 

The background evolution of dark energy can be parameterized by its equation of state 
wjc = px/px- The evolution of the energy density of dark energy is given by solving 



where the prime denotes the derivative with respect to the conformal time and TC = a' /a 
is the conformal Hubble parameter with a being the scale factor. As far as the background 
evolution is concerned, this equation is enough to specify the effect of dark energy on the 
evolution of the universe. 

However, as mentioned in the introduction, dark energy component itself can fluctuate 
and affect cosmic density fluctuations. For the fluctuation of dark energy fluid, we have 
to specify more in addition to the equation of state. Here we briefly summarize what we 
need to consider a general dark energy. In the rest of the paper, we follow the notation 
of Ref. [13] and work in the synchronous gauge whose metric perturbations are denoted 
as h and rj unless otherwise stated. We also sometimes make use of the gauge-invariant 
gravitational potentials, which appear as metric perturbations in the conformal Newtonian 
gauge, to discuss the effects of fluctuation of dark energy. 

To specify the nature of dark energy perturbation, we need to give its pressure pertur- 
bation and anisotropic stress. For the pressure perturbation, it may be useful to separate 
it to adiabatic and non-adiabatic parts. When the dark energy fluid is adiabatic, the evo- 
lutions of pressure perturbation can be specified with the adiabatic sound speed c 2 a which 
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is written as 

Then pressure perturbation is given as 5p = c 2 5p. If wx < and the time evolution of wx 
is small compared with the Hubble expansion rate, c 2 can be negative and the adiabatic 
pressure can not support density fluctuations. However, in general, non-adiabatic pressure 
fluctuation may arise. Such a degrees of freedom can be parameterized with the so-called 
(effective) speed of sound which is defined as 

2 _ S Px 



(3) 

rest 



Spx 

This quantity is usually specified in the rest frame of dark energy and we follow the 
convention here. A famous example of dark energy models with non-adiabatic pressure 
fluctuation is a scalar field, called quintessence, which has the speed of sound (? s — 1. Other 
models such as k-essence, a scalar field model which can have a non-canonical kinetic term, 
can have different values for the speed of sound. In the adiabatic case, c 2 and c 2 a coincide. 

In fact, to consider a general fluid model for dark energy, the speed of sound is not 
enough to specify its nature. We still have to determine an anisotropic stress for dark 
energy. Including the anisotropic stress, perturbation equations for density and velocity 
fluctuations are 



S'x = -(1 + w x ) [k 2 + 9H 2 (c 2 s - c 2 a )] ^ - m{c 2 s - w x )5x - (1 + w x )^ (4) 



0' x = -Ui\ - 3c 2 )6 x + t^-6 x - k 2 a x (5) 

1 + w x 



where 5x,0x and ax represent density, velocity and anisotropic stress perturbationa #1 l 
for a general fluid, respectively. The anisotropic stress can be presented as the viscosity 
and damps velocity perturbation in sheer-free frame. Thus, in addition to the sound 
speed of dark energy c 2 , there is another freedom to specify the anisotropic stress ax of 
dark energy. In fact, it may be possible to choose a x demanding that it is consistent 
with viscous damping of velocity perturbation in sheer-free frame. However, in this paper 
we follow a phenomenological approach of Ref . [TQl [12] . Considering the transformation 
property between frames, ax should be gauge invariant thus its evolution can be given by 
the gauge-invariant combinations of some metric and velocity shear. Taking into account 
the dissipation, the anisotropic stress can be obtained by solving [101 [12] 



a > +3n A ax = *(e x + K + 3v >) . ft;, 

Wx o \ 2 



# x The relation between ax here and ir g in Ref. [10] is 

2 w x 
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Here a is the viscosity parameter which specifies the nature of an anisotropic stress for 
dark energy. In fact, in some literatures, another parameter c 2 is is called the viscosity 
parameter, which is related to a as 

c 2 

a = ^xi^. (7) 
l + w x 

However, in this paper, we use a as the viscosity parameter. It should be mentioned that 
a must be kept positive to preserve the physical effects where the viscosity damps the 
perturbations. Thus, in the following we consider the case with a > 0. 

We modified CAMB code [H] to include a general dark energy fluid to calculate CMB 
and matter power spectra. We also checked our results by modifying CMBFAST [15] and 
found that they both give the same results with good accuracy. 



3 Effects of dark energy fluctuation 

In this section, we discuss the effects of dark energy fluctuation, paying particular atten- 
tions to the anisotropic stress of dark energy. Since the dark energy can be a dominant 
component of the universe only at late times, the effects of dark energy fluctuation is 
significant on large scalescj In particular, as we will show in the following, the effects of 
the anisotropic stress are significant at low multipoles of CMB power spectrum, through 
the late time integrated Sachs- Wolfe (ISW) effect. Since the ISW effect is driven by the 
change of the gravitational potential, we start our discussion with investigating the gravi- 
tational potential by changing the viscosity parameter a. In fact, as it has been discussed 
in the literature [TBI HTT]. the non-adiabatic pressure fluctuation, which is usually charac- 
terized by the effective speed of sound c 2 , can also affect the late time ISW effect. Thus 
we also discuss the effects of the non-adiabatic pressure fluctuation along with that of the 
anisotropic stress. 

Since we are going to discuss the effects on the ISW effect, first we show time derivative 
of the gravitational potential at present time in Fig. [Has a function of wave number k 

for several values of a and c 2 . Here and hereafter gravitational potentials are normalized at 
super horizon scales such that ^ = —10/ (4i? y + 15) with R u = p v j (p 7 + p u ) being neutrino 
fraction, and the conformal time coordinate is in units of Mpc following the convention. 
The equation of state for dark energy is assumed as wx = —0.8 in this figure. When we 
vary a, c 2 is assumed to be and vice versa. For illustrative purpose, other cosmo logical 
parameters are fixed as Q m h 2 = 0.128 (dark matter density), Q^h 2 = 0.022 (baryon 
density), h = 0.73 (hubble parameter), r = 0.09 (optical depth to the last scattering 
surface) and n s = 0.96 (spectral index of primordial density fluctuation), which are mean 



# 2 The background evolution of dark energy can affect the expansion of the universe, thus the peak 
positions of CMB power spectrum can depend on the equation of state of dark energy. However, as far as 
the fluctuation of dark energy concerned, its effects can be large only on large scales unless the equation 
of state dark energy is close to at early times. 
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values for ACDM model from the analysis using WMAP3 data alone [18] . In this paper, we 
assume that the universe is flat and no running of scalar spectral index and no tensor mode 
contribution. Unless otherwise stated, we use the above cosmological parameters except 
for those of dark energy. As seen from the figure, the change of the gravitational potential 
increases as a and c 2 becomes larger. Since a represents the strength of the anisotropic 
stress, which damps velocity perturbations, the increase of a drives the change of the 
gravitational potential. Thus we obtain larger ISW effect as a increases. As for the effects 
of the non-adiabatic pressure fluctuation, c 2 determines the sound horizon s = J c s dt/a 
under which dark energy can be assumed as "smooth" component. When dark energy is 
considered to be smooth, fluctuation of dark energy can act as stabilizing force for the 
gravitational potential, which leads to the decay of the gravitational potential. Above the 
sound horizon, dark energy can cluster, thus only gravitational force can be relevant to 
the evolution of the gravitational potential. The transition between "smooth" regime to 
clustered regime occurs at larger scales as c 2 increases. Thus the larger values of c 2 cause 
the more late time ISW effect, which can be observed in the figure. 
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Figure 1: The value of $' + today is shown as a function of k for several values of a 
(left) and c 2 s (right). The equation of state is fixed as wx = —0.8 in this figure. 



For the effects of the non-adiabatic pressure fluctuation, the behavior of the gravi- 
tational potential differs depending on c 2 only around the transition regime but not at 
smooth or clustered limit. However, for the anisotropic stress, the ISW effect depends on 
a on large scales. To see this clearly, we show the present value of $ + \I> as a function 
of k in Fig. [2j For the case with a being fixed and varying c 2 , it can be clearly seen that 
the behaviors on large/small scales are the same and they differ only at transition regime. 
However, for the case with c 2 being fixed and varying a, corresponding lines depend on the 
value of a even on large scale limit. Since the parameter a determines the magnitude of 
the anisotropic stress ax through Eq. flS]), the change of the gravitational potential is also 
affected by the magnitude of a. On the other hand, c 2 only determines the sound horizon 
which differentiates clustering property of dark energy, thus the magnitude of c 2 itself is 
not directly relevant to the amount of the change of the gravitational potential. This is 
one of the differences between the effects of c 2 and a although they affect fluctuation on 
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large scale in a similar way. 

Furthermore, it is well-known that $ and \& coincide when there is no anisotropic 
stress, in other words, $ — \I> is sensitive to the anisotropic stress perturbation a. This 
can be understood from the perturbed Einstein equation, as 

k 2 ($-y) = 12irG N a 2 (p + p)o- (8) 

where the quantity in the RHS is for the total matter. In the left panel of Fig. [3l the 
values of $ — \I/ at the present time are shown for several values of a as a function of wave 
number. As seen from the figure, |$ — \P| becomes larger as a increases on large scales, 
which can be easily understood since larger values of a means large anisotropic stresses 
from dark energy. Another point which should be noticed is that the effect of it becomes 
larger on large scales. For the generalized dark energy with non-zero a, the viscosity acts 
against the gravitational collapse once the modes come across the horizon (right panel of 
Fig. [3]). Therefore the perturbations including the anisotropic stress a have damped more 
at smaller scales, which makes the effect significant only at large scales at present. 

Now we are going to discuss the effects on CMB TT power spectra. In Figs. H] and 
El the CMB TT power spectra are depicted for several values of a and c 2 s . The value 
of the equation of state wx are found in the captions. As noted above, since the dark 
energy can be dominant component of the universe only at late times, its fluctuation 
affects the CMB power spectra on large scales. Thus we only show low multipole region in 
the figures. As seen from Fig. HJ Ci increases as a increases for the case with wx = —0.8 
and (? s = 0. This is expected from the argument on the ISW effect above. However, 
when the equation of state is less than —1, the tendency is the opposite, i.e., as the a 
increases, the large scale power decreases. This can be understood by looking at the source 
term in the perturbation equation. The anisotropic stress can affect the evolution of the 
gravitational potential through the change of density and velocity perturbations. Such an 
effect is induced by the source term in Eq. fll]) which depends on the prefactor 1 + wx- 
Thus the effect should be different depending on the sign of the factor 1 +wx- Hence the 
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Figure 3: (Left) Scale dependence of $ — ^ at the present time are plotted assuming 
several values for a. In this figure, c 2 are set to be 0. (Right) Time evolutions of ax for 
viscosity parameters a = 5 (black lines) and a = 1 (blue lines). Different lines correspond 
to the evolutions of different Fourier modes as indicated in the figure. The larger viscosity 
parameter makes the evolution of a after horizon crossing the slower during the matter 
dominated era. After the dark energy dominates the energy density of the universe, the 
perturbation starts to be erased. 



tendency becomes the opposite for the cases with wx > — 1 and < — 1. 

When the value of the sound speed is fixed as c 2 — 1, the effects of the anisotropic stress 
become insignificant for any value of wx- For the case with c 2 — 1, most scales relevant 
to low multipole region where dark energy fluctuation can be important are in smooth 
regime. As shown above, the anisotropic stress can affect the gravitational potential, at 
least superficially, similarly to smoothed dark energy. Hence when dark energy can be 
assumed as a smooth component, the effect of the anisotropic stress becomes irrelevant 
since dark energy is already smoothed. 

For comparison, we also show the same plot for the effects of the non-adiabatic pressure 
fluctuation. In Fig. [5j CMB TT power spectra are shown for several values of c 2 s fixing 
a and other parameters. Since the effects of the anisotropic stress and non-adiabatic 
pressure fluctuation are similar in the ISW effect, the tendency should be also the same 
as that of the effects of a. As anticipated, in Fig. [SJ we can see almost the same behavior 
of low multipole region for the case with c 2 being varied, as for the case with a being 
varied. Thus we can expect that there are some degeneracy between the effects of a and 
c 2 on the fluctuation at low multipoles of CMB. As discussed above, they affect the ISW 
effect which enhances/reduces the power at low multipoles depending on the strength of 
these effects characterized by a and c 2 . Although the origin of these effects are different 
in nature, they seem to affect low multipoles of CMB in the same manner. To see this 
quantitatively, we show the contours of constant ratio of C\ = 1(1 + 1)C; at between I = 2 
and the first peak in Fig. [6j In the figure, the equation of state of dark energy is assumed 
as wx = —0.8. Other cosmological parameters are also fixed as the values already noted 
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above. Notice that higher multipole region I > (9(100) is unchanged by varying the values 
of c 2 and a. The ratio represents how low multipoles are enhanced/reduced by the effects 
of the anisotropic stress and non-adiabatic pressure fluctuation. The figure indicates that 
some combination of c 2 and a give the same effect on large scales. Thus when we study 
the constraint from observations, we can expect some degeneracy between c 2 and a. We 
discuss this issue in the next section. 
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Figure 4: CMB TT power spectra for several values of a for the cases with wx = —0.8, c 2 = 
(top left), wx = — 1-2, ci? = (top right), wx = —0.8,c 2 = 1 (bottom right) and 
Wx = —1-2, c 2 = 1 (bottom right). 



Here we should make a comment on a nontrivial cancellation between the ISW effect 
and the cross correlation between the ISW and Sachs- Wolfe (SW) effect. Although we 
argued that larger values of a give the larger ISW effect for the case with wx > — 1, the 
lines of a = 1 and 5 in the top left panel of Fig. |4]do not obey this rule. In fact, this is 
because the ISW and SW effects can nontrivially cancel those effects in some cases. We 
discuss this issue briefly. For this purpose, we make use of the transfer function on large 
scales. Formally, C\ is written as 

/dk 
—P n {k)\^{k)\ 2 . (9) 

where Pn{k) is the primordial power spectrum and Ai(k) represents the transfer function 



8 



3000 p 

2500 
2000 
1500 
1000 

500 - 

- 
1 



=0.001 




10 

multipole / 




10 

multipole / 



■3 



3000 
2500 
2000 
1500 
1000 
500 




-=o.ooi 

2 




10 

multipole / 




10 

multipole / 



Figure 5: CMB TT power spectra for several values of c 2 . for the cases with wx = —0.8, a 



(top left), wx = —1.2, a = (top right), wx = —0.8, a 
wx = —1.2, a = 1 (bottom left). 



1 (bottom right) and 



for a multipole moment I at present time. On large scales, Ai(k) are mostly determined 
by the contribution from the SW and ISW effect, thus it can be written as 

A t {k) = Af w (k) + A] sw (k) . (10) 

In particular, the contribution from the ISW effect is given by 

A] sw (k) = J d V e-^ + V) Mk( v - m )] , (11) 

where r is the optical depth along the line of sight, $ is the gravitational potential and 
ji is the spherical Bessel function. We plot |AJ, SW | 2 as a function of wave number k for 
several values of a in Fig. [7| for the case with wx = —0.8. As seen from the figure, | A2 SW | 2 
becomes larger as a increases. In fact, as mentioned above, the quadrupole Ci for a = 1 is 
larger than that for the case with a = 5 although the ISW effect is larger for a = 5. Thus, 
at first sight, it looks contradictory. However, this comes from a non-trivial cancellation 
between the contributions from the ISW effect | Z\^ SW 1 2 an d the cross correlation between 
the SW and ISW effects A^ w x Af w . In Fig. El we also plot A^ w x A ! 2 SW as the same 
manner as Fig. [7J When a is larger, the cross correlation becomes also large in amplitude 
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Figure 6: Contours of the ratio C 2 /Ci st are plotted in the c 2 s — a plane for the case with 
Wx = -0.8. 

as seen from the figure. However, notice that, since the cross correlation of the SW with 
the ISW effects can give a negative contributions to the transfer function. Thus, for some 
cases, it gives smaller quadrupole amplitude as a whole effect. This is the reason why 
the low multipoles for the case with a = 1 is larger than that for a = 5 even though 
the contribution from the ISW effect monotonically becomes larger as increasing a. For 
reference, we also show the same plot for the case with wx = —1.2 on the right panel in 
Figs. [3 and El 

4 Constraints from current cosmological observations 

In this section, we show the constraints from current cosmological observations on the 
energy density and equation of state of dark energy . 

To include the effects of non-adiabatic pressure fluctuation and anisotropic stress of 
dark energy, we calculate theoretical angular power spectrum by CAMB code modified ac- 
cording to the details in section 2. Because these effects change the power at larger angular 
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Figure 7: Contribution from the ISW effect on the quadrupole [A2 SW ] 2 are shown for 
several values of a for the cases with wx = —0.8 (left) and —1.2 (right). 
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Figure 8: Contribution from the ISW effect on the quadrupole A 2 SW x A^ w are shown for 
several values of a for the cases with wx = —0.8 (left) and —1.2 (right). 



scales, the amplitude and the shape of the power spectrum at low multipoles can be used 
in principle to put constraints on these fluctuation properties. In fact, however, we cannot 
find any significant limits on them because of the large errors due to the cosmic variance at 
low multipoles. Therefore, we conclude that the current cosmological observations cannot 
put tight constraints on the a or (? s parameters. (The forecast of observational constraint 
on c 2 s from planning galaxy surveys can be found in [T9].) 

Even though one cannot give limits on these parameters themselves, it does not mean 
that they are irrelevant when considering the constraints on the nature of dark energy. 
It can happen that one obtains different limits on the cosmological parameters such as 
equation of state of dark energy wx if one takes different parameters on the fluctuation 
property of dark energy such as effective sound speed c 2 s and/or anisotropic stress a. Many 
observational proposals aim to determine the equation of state parameter of dark energy 
w x as precise as possible, even up to ~ 1% level [20], because the information about wx 
is important to pin down a model for dark energy, and also to determine the future of 
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universe [21 J. It is therefore important to investigate how much different assumptions on 
the fluctuation properties of dark energy can lead the different observational limits on the 
equation of state of dark energy. 

To study this explicitly, we made likelihood analysis to put constraint on wx with 
several different assumptions on c 2 and a. The likelihood functions we calculate are based 
on WMAP three year data [221 6] and the 2dF Galaxy Redshift Survey [21]. To make 
analysis general enough, one has to vary the other cosmological parameters which also 
affect the shape of the CMB and matter power spectra. In order to do so effectively, 
we follow Markov chain Monte Carlo (MCMC) approach [25J, and explore the likelihood 
in seven dimensional parameter space, namely, 9 (the ratio of the sound horizon to the 
angular diameter distance at last scattering), A s (amplitude of primordial perturbation), 
n s , r, and w X - Here the Hubble parameter h is derived from 9 [26J and the 
dark energy density such that the universe is spatially flat. 

In the left panel of Fig. [HI marginalized 1 and 2cr allowed regions in the Q m — w x 
plane from WMAP three year data are shown for the case with c 2 = and a = 0. For 
comparison, the case with c 2 = 1 and a = is also shown in the figure. As seen from the 
figure, the constraint on w x for the case with c 2 = and a = becomes severer compared 
to the other case. This is due to the fact that the ISW effect becomes significant as wx 
decreases to negative value when wx < — 1, which makes the fit to the data worse at 
low multipole region. This results indicates that the perturbation nature of dark energy 
fluid affects the determination of the equation of state in some cases. The marginalized 
one-dimensional probability distribution of wx is also shown in the right panel of Fig. [91 
We found that the confidence region of wx can differ by ^ 10% between the models with 
(a,c 2 s ) = (0,0) and (0,1). 

In the left panel of Fig. \TU\ we show marginalized 1 and 2a allowed regions in the 
Q m — wx plane obtained by combining the information from the matter power spectra by 
2dF galaxy survey with WMAP three year data. Although the data reduce the allowed 
region in the Q m — wx plane, it can still be found that there is a significant difference 
between the models with different assumptions on the fluctuation properties of dark energy. 
Similar analysis and conclusion are done and derived by [33] for a fluid dark energy model 
parameterized by w and c 2 , by cross-correlating the map from WMAP three year data 
with that from NRAO VLA sky survey data (see also [17J). The present work differs from 
them in that we analyzed different observational data set, included viscosity parameter a, 
and relaxed the other cosmological parameters to be varied while they fixed cosmological 
parameters except for the dark energy parameters. 

In Fig. [HI marginalized 1 and 2a allowed regions are shown for the case with (c 2 , a) = 
(0,0.6) (left) and (1,0.6) (right). In fact, large value of c 2 /a makes any dependence of 
Wx at low multipoles insignificant. Furthermore, we found that when a ~ 1 the result 
is insensitive to the value of c 2 because of the fact that the fluctuation amplitude at low 
multipoles is independent from c 2 in such cases (see Fig. [6]). Thus in this case, the 
constraint on Q m — wx plane is almost unchanged. We also show the same plot for the 
case where we use the data from 2dF in addition to WMAP three year data in Fig. [121 As 
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Figure 9: (Left) Allowed regions of la (red region) and 2a (purple region) are shown 
for the case with c 2 s = and a = 0. For comparison, allowed regions for the case with 
Cg = 1 and a = are also shown in black dashed lines. (Right) Marginalized probability 
distribution of the equation of state parameter wx with different assumptions on the 
fluctuation properties of dark energy. Black line corresponds to the case with (c 2 sl a) = 
(1, 0) and red dashed line with (c 2 , a) = (0, 0). 
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seen from the figure and easily expected, the constraints become severer but the differences 
among different assumptions on c 2 s and a are small. 




Figure 11: Allowed regions of la (blue region) and 2a (cyan region) are shown for the case 
with cj. = and a = 0.6 (left) and cj. = 1 and a = 0.6 (right). For comparison, allowed 
regions for the case with c 2 s = 1 and a = are also shown in black dashed lines. 




5 Conclusions and discussions 

We discussed the effects of a generalized dark energy on cosmic density fluctuations and 
studied cosmological constraints on dark energy in that setting. When one considers 
a general type of dark energy, one can also include a possible anisotropic stress of dark 
energy, which has not been investigated much in the literature. We studied the effects of the 
anisotropic stress characterized by the viscosity parameter a along with the non-adiabatic 
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pressure fluctuation usually described by the speed of sound c 2 . First we discussed the 
effects of them on the gravitational potential and argued how the values of a and c 2 affect 
it. As shown, in section [HJ the anisotropic stress of dark energy can affect the ISW effect in 
a similar way with the non-adiabatic pressure fluctuation. We also explicitly showed that 
there is some degeneracy between the effects of a and c 2 on the power at low multipoles 
of CMB angular power spectrum. 

We also explore constraints on dark energy from recent cosmological observations as- 
suming a generalized dark energy. As shown in Fig. [91 the constraint on the equation 
of state can be changed depending on the parameters c 2 and a which describes the per- 
turbation nature of dark energy. This is because fluctuation on large scales are affected 
by the nature of dark energy fluctuation, i.e., depending on a and c 2 , through the ISW 
effect. This analysis may indicate that cosmological constraint on dark energy such as 
the equation of state should be carefully done by taking the perturbation nature of dark 
energy into account. 

Acknowledgments: KI acknowledges the support by a Grant-in-Aid for the Japan So- 
ciety for the Promotion of Science. 
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